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We develop a minimal multiorbital tight-binding model with realistic hopping parameters. The model breaks 
the symmetry of the tetragonal point group by lowering it from C4 to Z?2d, which accurately describes the 
Fermi surface evolution of the electron-doped BaFe2- ;E Co a ;As2 and hole-doped Bai_ H K y Fe2As2 compounds. 
An investigation of the phase diagram with a mean-field t-U -V Bogoliubov-de Gennes Hamiltonian results in 
agreement with the experimentally observed electron- and hole-doped phase diagram with only one set of t, U 
and V parameters. Additionally, the self-consistently calculated superconducting order parameter exhibits s - 
wave pairing symmetry with a small d-wave pairing admixture in the entire doping range, which is the subtle 
result of the weakly broken symmetry and competing interactions in the multiorbital mean-field Hamiltonian. 

PACS numbers: 78.70.Dm, 71.10.Fd, 71.10.-w, 71.15.Qe 



Introduction. — The discovery of the iron-based supercon- 
ductors attracted much experimental and theoretical attention 
leading to many microscopic model studies ifTT- tlOl . Of the 
various effective microscopic models, the one proposed by 
Zhang [ 1 ] accounted for the effects of the upper and lower an- 
ion atoms (above and below the Fe plane). Zhang's minimal 
model successfully described the behavior of the collinear an- 
tiferromagnetism (C-AFM) ifTTT - TPIl and the competition with 
superconducting (SC) order in the electron-doped (e-doped) 
part of the phase diagram when substituting Fe with Co atoms. 
However, the model failed to give account of the hole-doped 
(/i-doped) part of the phase diagram when doping on the Ba 
site. Very recently, Hu and co-workers (7] proposed a Hamil- 
tonian with S4 symmetry to clarify the local symmetry break- 
ing of the underlying electronic structure through orbital or- 
dering in the Fe-122 family. Both models by Zhang [JT| and 
Hu share a very important ingredient, that is, the breaking 
of the C4 symmetry. So far both models have been applied 
only to the e-doped region of the phase diagram 12 El- Al- 
though both models give a qualitative picture of competing 
order in the system, caution must be taken when the evolu- 
tion of the low-energy quasiparticle states and Fermi surface 
topology are considered. These can impose even stricter con- 
straints on a given model: (i) e-doping: The disappearance of 
the hole pockets at the T point of the Brillouin zone (BZ) is the 
key feature of e-doped compounds BaFe2- K Co x As2 lfl4ifi"8l . 
(ii) /i-doping: The nested large electron pockets around the 
M point of the BZ evolve into a set of four small clover-like 
hole pockets for the end member KFe 2 As 2 , which is identi- 
fied as a Dirac cone |T9]42TI . A very similar feature is also 
found in several density functional theory (DFT) calculations 
E] |9] |22l ES. Note that this is an intrinsic feature of the 
bare band structure and is not related to emergent spin-density 
wave order I1TT1I241 . 

Aided by experiments and DFT calculations, we develop 
a minimal tight-binding model with improved normal-state 
band structure parametrization to account for conditions (i) 
and (ii) . The crystal structure of BaFe2As2 (ThCr2Si 2 type 
structure) and its electronic structure are now well understood 




FIG. 1. (color online) (a) Schematic picture of the Fe d xz and d yz or- 
bital overlap through the As atoms. In panel (a), the red (blue) lobes 
represent the d xz (d yz ) orbitals. The white (black) spheres represent 
the Fe atoms on the A (B) sublattice and the yellow spheres represent 
the upper (lower) As atoms. This configuration shows a 2x2 (4-Fe) 
unit cell within the solid thick lines. The shaded region is the 2-Fe 
unit cell. The 2-Fe unit cell is used in panel (b) to construct the tight- 
binding model, where hopping through the As atoms is ignored. In 
panel (b), the solid black line stands for nearest-neighbor (INN) intra 
(ti) and inter (ts) hopping parameters, the solid red (blue) line stands 
for 2NN hopping t2 (U), the dashed red (blue) line is for second- 
nearest-neighbor (2NN) hopping t3 (t4), and the dotted black line 
stands for the third-nearest-neighbor (3NN) hopping te. 

and the relevant orbitals for atomic bonding have been identi- 
fied E5ll . Within a quantum-chemical framework of bonding 
one realizes that the upper (lower) anion atom gives rise to 
different overlap between the Fe 3d xz (3d yz ) orbital with the 
4p orbitals of the As atom. A schematic picture of the overlap 
between these orbitals is shown in Fig.[TJa). The challenge in 
the study of superconductivity in the Fe -based 122 family is 
how to simultaneously satisfy the observations (i) and (ii). 

In this Letter, we show that in order to model the phase 
diagram, one needs to introduce three key conditions critical 
for the second-nearest neighbor (2NN) hopping terms. I: On 
the same orbital, we break C4 symmetry between the A and B 
sublattices of the Fe atoms (t 2 7^ £3), see Fig.[TJb). II: In the 
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FIG. 2. (color online) The electronic band structure of model Hamil- 
tonians Hi (a) and H2 (b) along high-symmetry directions in the BZ. 
The insets show the significant difference between and ¥l 2 near 
the M point. The horizontal line in (a) and (b) indicates the Fermi 
level Ef for half filling (n = 2.0). Panel (c) shows the Fermi sur- 
faces for model for n = 2.23 (blue, fi=-0.76), n — 2.0 (green, 
^=-1.244), n = 1.8 (cyan, /i=-1.65), n = 1.5 (red, /i=-2.25). Panel 
(d) shows the FSs of model H 2 for different doping levels n = 2.38 
(blue, ^=-0.055), n = 2.0 (green, ^=-0.89), n = 1.8 (cyan, fi=- 
1.25), n = 1.5 (red, ^=-1-745). 



same sublattice, we break the degeneracy between the d xz and 
d yz orbitals by adding a twist. Ill: We include the 2NN inter- 
orbital hopping term £4. The combined effects of conditions I 
and II reduce the symmetry to D^d- For a detailed explanation 
see the Supplemental Material (SM) l26l . 

Theory.— The original set of four hopping parameters, 
ti_4, in Zhang's model accounted only for conditions I and 
III. In our new model, we include the symmetry condi- 
tion II and thus extend the set to six hopping parameters 
ti_6' For the calculations we choose the 2-Fe unit cell, 
as described in Fig. fTTa), to construct the basis function ip 
= (cai, CA2, cbi, cb2) and its Fourier transform ciua = 
^= ^ k CMa(k) x exp (iR,;Af • k) for the normal state k- 



x exp 

space. The kinetic energy is then given by the minimal hop- 
ping Hamiltonian H l — ^ k V^(k) ip(k) in the 2-Fe unit 
cell of the 1 st BZ, -vr < k x (k y ) < n, where M=(A,B) and 
RiM is the position of the Fe atoms in the A (B) sublattice, 
a=(l,2) stands for d xz I d yz orbitals, and 
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which satisfies conditions I through III 12711 : 

Hi ; iAi = iA2 = i , £bi = iB2 = i v , 

with tx- 6 = (-1, -0.4, 2, -0.04, 0, 0). 

€ai = iB2 = i H , iA2 = iBi = i v , 

with t 1 _ 6 = (-1, 0.08, 1.35, -0.12, 0.09, 0.25). 



(2) 



First, we compare Zhang's minimal model H|, which sat- 
isfies only conditions I and III, with our extended model H|, 



Here, we defined the dispersion functions: £ H = 
2t2COs(k x ) + 2t 3 cos(fc a ) + 4t 6 cos(k x ) cos(ky), £ v = 
2t 3 cos(k x ) + 2t2COs(k y ) + 4t 6 cos(k x ) cos(k y ), ^i 2 = 
2t4COs(k x ) + 2t4COs(k y ), £ t — 4ii cos(k x /2) co$(k y /2), 
^ c = At§ cos{k x /2) cos(k y /2). Fig.[TJb) shows the schematic 
picture of H^. The key difference between both models is the 
orientation of each orbital, which is identical in the A and B 
sublattices in Zhang's model, while it has a 90° relative ro- 
tation (twist) in the new model. The parameters for H\ are 
taken from Ref. 2 In the case of H\ we determine the hop- 
ping parameters by comparing the calculated FS topologies to 
the angle-resolved photoemission (ARPES) experiments |[T9l 
for electron fill factors n = 2.0 (half filling), n = 1.8 (opti- 
mal /i-doping) and n = 1.5 (KFe2As2). In Fig. [2ja) and (b) 
we plot the band dispersion for H* and H 2 , respectively. Note 
the significant difference between the electronic dispersions of 
these Hamiltonians, which is seen in the insets of Fig. [2| a) and 
(b), where we enlarge the relevant area around the M point in 
the BZ. In the case of (b) a linear dispersion (Dirac cone) can 
be found for model H|, while none exists for model in (a). 
As will become clearer, this key difference between Zhang's 
H\ model and the new H 2 model is a direct consequence of 
condition II. 

In the next step, we construct the minimal mean-field 
Hamiltonian in the weak-coupling regime. Since the C-AFM 
will enlarge the real-space unit cell, we need to choose the 
4-Fe unit cell configuration as shown in Fig. JTJa), with the 
real-space Hamiltonian, H = H* + H A + H mt . Where H* 
is the hopping term, H A is the 2NN intra-orbital pairing in- 
teraction and H" lt is the interaction term, which includes the 
Coulomb interaction U and Hund's coupling Jh ■ For details 
of the construction of the Hamiltonian see the SM |26|. 

Model results.— We calculate the phase diagram for both 
Hamiltonians H % =H\ + H lnt + H A with i = 1,2 to in- 
vestigate the stability of the reported C-AFM, SC, and co- 
existing C-AFM/SC phases as a function of doping param- 
eter. The evolution with doping is described by the same 
set of hopping, interaction, and pairing parameters for each 
model across the entire doping range. For model Hi we use 
(U, J H , V) = (3.4, 1.3, 1.2), while for H 2 the set of parame- 
ters (U, J H , V) = (3.2, 0.6, 1.05) is used JH. 

Figure |3ja) shows the C-AFM and SC order parameters at 
nearly zero temperature (T = 10~ 4 ) for models Hi (top) 
and H2 (bottom). We confirm that the phase diagram of e- 
doped compounds is equally well described by both Hamil- 
tonians in agreement with experiments |fl3HT8l l29l . How- 
ever, the situation is markedly different on the /i-doped side. 
Here only H2 is capable of describing experiments ll20l |2TI 



3 



0.14 

0.12 

0.10 

0.08 

0.06 

0.04 

0.02 

0.00 
0.14 

0.12 

0.10 

0.08 

0.06 

0.04 

0.02 

0.0^ 



(a) 



|;n,| 
~ A. 



...B 


2 










*~ ~* \ m i\ 






!i" 

i 


1 










f 


» - 

■■<!•- 










i 




1 














i 




t 














i 
j 













5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 

n 




(c) 



:A rf /A 



1.5 1.6 1.7 1.8 1.9 



FIG. 3. (color online) The phase diagrams of the C-AFM (m.) 
and SC (A s ) order parameters of models Hi (top) and H2 (bottom) 
at zero temperature (T = 10 -4 ) are shown in panel (a). Panel (b): 
Temperature dependence of the calculated order parameters of model 
H2. The suppression of rru is visible in the coexistence region with 
A s at filling n = 1.9. Panel (c): The ratio of the SC order parameters 
Ad/ A s for 2NN pairing is shown for Hi. The d-wave admixture is 
of order 8% over the entire ft-doping regime (n = 1.5 — 2.0). 



by correctly accounting for the e-h asymmetry and the exis- 
tence of a strong SC phase at low electron filling (1.6 < n < 
1.85). In Fig. |3jb) the finite temperature self-consistent cal- 
culations of the C-AFM and SC order parameters are shown 
for H 2 at /i-doping values n = 1.5, 1.6, 1.7, 1.8, 1.9 and 
2.0 (half filling). We also find that in the coexistence C- 
AFM/SC phase both orders compete for phase space. The 
competition leads to a marked suppression of rrii, when A s 
nucleates at a lower temperature, see Fig. |3jb) for fill fac- 
tor n = 1.9. From the self-consistent calculations of the 
phase diagram, we extract the maximum gaps: C-AFM gap 
rrn/T N w 0.116/0.211 = 0.549 at n = 2.0, and the SC 
gaps for e-doping Af p /T c « 0.11/0.049 = 2.23 at n = 2.1, 
and /i-doping AfP/T c w 0.14/0.066 = 2.16 at n = 1.8 
and AfP/T c « 0.033/0.015 = 2.21 at n = 1.5. Here we 
introduced A e s xp = 4A S as the experimentally determined 
tunneling gap. The SC gap ratios are in reasonable agreement 
with experimental reports for various doping values, ranging 
from values associated with weak- to strong-coupling pairing, 
1.5 < A ex P/T c < 3.8 I30H331 . Fig. gc) shows the ratio 
between Ad and A s , which is roughly 8% over the entire h- 
doping region. Although, the s gap symmetry is the most 
widely accepted pairing symmetry for Fe-122 based SCs, we 



always find a small admixture of c?-wave symmetry in our self- 
consistent mean-field calculations. A weak admixture of the 
d-wave channel is not unexpected, since H 2 breaks the C4 
symmetry of the point group. For a lattice model with D 2 d 
symmetry both s- and d-wave belong to the same representa- 
tion and are allowed to mix. This admixture is also found in 
Zhang's model, H\, however, there the C4 symmetry break- 
ing is only due to condition I, while for model H 2 it is related 
to condition I and the 2NN hopping terms of condition II. 

The pairing symmetry can be discussed more systemat- 
ically by how it affects the spectral function observed in 
ARPES experiments. Since we are interested in the en- 
tire doping range of the phase diagram, we focus only on 
the H 2 model. The discussion can be further simplified 
by neglecting the coexistence region of C-AFM/SC. In that 
case, we can downfold the BdG Hamiltonian from the 4-Fe 
unit cell onto the 2-Fe unit cell configuration, see Fig. [TJa). 
The order parameters A s and A^ were previously defined 
on the lattice. Their Fourier transforms in the BZ of the 
2-Fe unit cell are A s (k) = 2A S \cos(k x ) + cos(fc y )] and 
Arf(k) = 2 Ad [cos(fca;) — cos(ky)]. We can further decom- 
pose them into sublattice (A,B) and orbital (1,2) contribu- 
tions: A A1 (k) = A B2 (k) = A s (k) + A d (k), A A2 (k) = 
Afli(k) = A s (k)-A d (k). 

Note these expressions are a consequence of the twofold 
D 2 d symmetry, thus giving rise to the s x 2 +y 2±d x 2 _ a 2-wave 
gap. The corresponding downfolded fc-space BdG Hamilto- 
nian attains the eigenvalues and eigenvectors for the spectral 
function of each orbital and sublattice, 



Aj a (k,w) 



n=l 



i< t (A0i^-£n+l< Q #)r^+£^ 



(3) 

where the total spectral function, measured in ARPES, is 
given by the sum A(\t, oS) = ^ Ai a (k, uj). An obvious 
question is what is the evolution of the spectral function with 
doping and how is it affected by the order parameters A s 
and Ad for filling factors n = 2.1, 1.8 and 1.5. In Fig. [4] 
the total spectral function and its partial weights are shown. 
The reduced D 2 d symmetry is obvious from plots of spec- 
tral weights Aai,Aa2,Abi and Ab2- By construction of 
the model, the 90° rotational symmetry breaking (twist) is re- 
vealed by the distinction between A^xiBi) an d Aj^mi) (right 
panels in Fig. |4j. Unfortunately, ARPES experiments cannot 
differentiate between them. However, one can compare the 
total spectral functions (left panels of Fig. |4j with ARPES ex- 
periments. The spectral weights are very small for n = 2.1 
and n = 1.8, indicating that quasiparticles are gapped. The 
hole pockets around the V point exhibit an isotropic gap for 
n = 2.1, while an anisotropic gap with fourfold modulation 
is found for n = 1.8. On the other hand, for n = 1.5 (cor- 
responding to KFe2As2), we observe eight quasiparticle hot 
spots at the corners of the large hole Fermi surface centered at 
r, as well as four small pockets near the M point. We believe 
that the hot spots indicate either the existence of nodal octet 
structure or highly anisotropic gap. 
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FIG. 4. (color online) Spectral functions at zero temperature at the 
Fermi level for fill factors n = 2.1 (a), n = 1.8 (b), and n — 1.5 
(c). In each panel the total weight of the spectral function is shown 
in the left column, while the right column shows the corresponding 
partial spectral functions for each orbital and sublattice. The total 
spectrum is the sum of all four partial spectra. Panel (c): The two 
arrows mark the peaks of the electronic hot spots responsible for the 
octet gap structure. The dashed line shows the nodes of A 3 (k) = 0. 



The pairing symmetry in KFe2As2 has been a recent topic 
of hot debate. The existence of only hole pockets has moti- 
vated earlier theoretical proposals of d-wave pairing symme- 
tries with gap nodes l34ll35ll . This would imply a change in 
the superconducting symmetry from s- to d-wave pairing as 
hole doping is increased. So far there is some experimental 
evidence for gap nodes [36-39]. However, more recent angle- 
resolved photoemission spectroscopy unveiled that KFe 2 As 2 
is a nodal s-wave superconductor with octet-line node struc- 
ture on the large hole pockets at k z = tt ll40l . Our model 
calculations are consistent with an octet nodal structure on the 
zone-centered hole pocket at k z = 0. However, the eight elec- 
tronic hot spots at the corners of the large hole Fermi surface 
in the spectral function do not rule out a highly anisotropic 
gap. Future ARPES measurements at k z = may resolve the 



current disagreement about the location of hot spots vs. nodal 
points between our results and experiments. 

Summary. — We have shown how the C4 symmetry break- 
ing, involving the As atoms below and above the Fe layer, 
leads to a natural extension of Zhang's minimal model for the 
Fe-122 superconductors. The new tight-binding model is in 
better agreement with experiments and density functional the- 
ory calculations over the entire doping range by incorporating 
an orbital twist. Within the weak-coupling theory of super- 
conductivity, the calculated phase diagram reproduces quali- 
tatively the e-h doping asymmetry reported in many experi- 
ments and the gap-over- T c ratios. In addition, we always find 
a small d-wave admixture of roughly 8% to a dominant s- 
wave SC order parameter. This admixture can give rise to 
quasiparticles in the spectral function resembling nodal ex- 
citations at electronic hot spots in the spectral function of 
KFe2As2- Finally, the new minimal model is computation- 
ally more efficient than similar five-band models for studying 
disorder effects in real space around impurities or magnetic 
vortices. 
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SUPPLEMENTAL MATERIAL 

Here we present the construction of the mean-field Hamil- 
tonian on a two-dimensional lattice employed in the main text. 

THE MEAN-FIELD BOGOLIUBOV-DE GENNES 
EQUATION 

For the mean-field Hamiltonian, H = E l + H A + H mt , 
the hopping Hamiltonian H f is expressed by 

H = (tiaju'Ci aa Cja'a + h.c) — [I ^""^ c\ aa Cj aa , (4) 

iajct'a iacr 

where i and j are site indices of each Fe site, a=(l,2) is the or- 
bital index, a = (t, I) is the spin index and \i is the chemical 
potential. The second term in H is the mean-field supercon- 
ducting (SC) pairing Hamiltonian, 

H A = J2 (A^cLaCW + h.c). (5) 

ijOLCr 

The third term in H is the on-site Coulomb interaction Hamil- 
tonian including Hund's coupling. Following Ref. 1 we write, 

+ (U' - J H ) 2_j ( n ia<y) n ia'<r, 

where n iaa = c\ aa c iaa and U'=U — 2 . After performing 
the Bogoliubov transformation of the electrons onto particle- 
hole excitations with c laa = u" QCT 7„ + a v™* a ^ n , the corre- 
sponding Bogoliubov-de Gennes (BdG) Hamiltonian can be 
written in matrix form and solved self-consistently, 

(7) 



THE D 2d INVARIANT SYMMETRY 

We now give the reason for our choice of the Did sym- 
metry of the point group describing the crystal structure of 
BaFe2As2, which is used for the construction of the kinetic 
Hamiltonian in Eq. (1). To begin with, we draw the three- 
dimensional (3D) structure of the basic Fe-As building blocks 
of BaFe2As2 in Fig. [5] The D^d symmetry is generated by the 
group elements C4 and ah as shown in Fig.[5Jb,c). Close in- 
spection reveals the point group symmetry D 2c i, because com- 



with particle-like (u" a ^) and hole-like (v" a ^) wave functions. 

The self-consistent mean-field equations for the SC order 
parameter and spin-up and spin-down occupations are given 
by 

n \ B / 

where 

(n iat )=J2K a t\ 2 f(En), (9) 

n 

and 

(n^) = ^|^| 2 [l -f(En)}. (10) 

n 

For simplicity, we consider second-nearest-neighbor (2NN) 
intra-orbital pairing only with pairing potential Vy=V. To fa- 
cilitate the discussion of physical observables and generating 
of the phase diagram, we define the staggered lattice magne- 
tization and the s-wave and d-wave projections of the order 
parameter on the lattice: 

m = j y^(( n ict) - {niaj}), (11) 

a 



A * = giV ^ A m+*><*' (12) 



Ad = g^H 6x6 y A v+«,«|- (13) 

i,8, a 

The neighbors of site i are reached by <5 = (x + y,x — y,— x + 
y, —x — y), with e x — S ■ x,e y — S ■ y; N is the number of 
Fe sites in the real-space lattice. We note that the x (y) axis is 
aligned with the short Fe-Fe bond direction. 

I 

bined fourfold rotation (C4) and mirror reflection (07,.) leave 
the crystal structure invariant. 

However, the D 2 d symmetry is obvious only for the 3D 
crystal structure, while the inclusion of the reflection opera- 
tion <Jh is not obvious for the two-dimensional (2D) model 
as shown by Fig. 1(a) in the main text. The question is how 
one can construct a D^d hopping Hamiltonian on a 2D lattice. 
The solution is as follows: The As atoms mediate the 2NN 
hopping (£3) through their p orbitals between the d orbitals 
of the Fe sites. As a consequence the upper (lower) As atoms 
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FIG. 5. (color online) (a) Schematic picture of D2d point group symmetry of the three-dimensional (3D) structure of the Fe-As building blocks. 
The colormap is: Fe atoms on the A (black) and B (white) sublattice; the As atoms above (green/yellow) and below (yellow/green), (b) Top 
view of panel (a); (c) application of symmetry operation C4 or oh on panel (b); (d) orbital ordering and tz hopping terms overlayed onto panel 
(b); (e) the resulting ah operation of panel (d). 



lead to effective hopping terms between the d xz (d yz ) orbitals, text, 
respectively. Finally, the 07, operation of the upper/lower As 

atoms can be mapped onto the exchange of the order of the 

d xz and d yz orbitals, see Fig. |5Jb). This corresponds to ex- 
changing the upper and lower panels of Fig. 1(b) in the main 
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